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Abstract 

The stable chromomagnetic vacuum for SU{2) Yang-Mills theory found 
earlier is shown to give a model for confinement in QCD, using Wilson loop 
and a linear potential (in the leading order) for quark- antiquark interaction. 
The coefficent k in this potential is found to be ~ 0.25 GeV^, in satisfactory 
agreement with non-relativistic potential model calculations for charmonium. 
At finite temperature, the real effective energy density found earlier is used 
to obtain estimates of the deconfining temperature agreeing reasonably with 
lattice study for SU{2). 



The economical definition of confinement of quarks in QCD is the 'area 
law' for the Wilson loop. The gauge invariant Wilson loop is 

W{C) = Tr P e-'^f '^^'^ ^^W^", (1) 

where P denotes the path ordering and t"' are the generators of the gauge 
group. We shall consider SU{2) Yang-Mills theory and choose the Savvidy 
[1] classical backgroound 

A^, = ; At = 5^'i-^,^,0). (2) 
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This choice solves the classical equation of motion JJ^^pt^^^ = Q, where 
D^^ = d^S"''^ + ge"'^'^A^^. The classical background corresponds to constant 
chromomagnetic field in the third color direction = H and this comes 
from the derivative terms in F'^^ = d^A'l — duA'^^ + ge°'^'^A^^Al. For this reason, 
the Savvidy ansatz (2) is called 'Abelian like'. So the classical background (2) 
is esssentially Abelian- like, taking values in the Cartan subgroup of SU{2). 

The use of Abehan-hke field strength can be understood from the idea of 
't Hooft [2] who proposed 'Abelian Projection'. This is a particular gauge 
fixing, breaking the gauge group SU{N) to the maximal torus subgroup 
H = f/(l)^~^. For SU{2), H = U{1). This is realized in a specific gauge 
called the 'Maximal Abelian Gauge'. In the continuum formulation, this has 
the form 

{d,6'^'' + ge^''Al)Al, = 0, (3) 

and the classical Sawidy background (2) satisfies (3). Numerical simulations 
on the lattice have found that the Abelian projected Wilson loop defined by 
A^ exhibits the 'area law' [3]. So (1) becomes 

^ (g-^f^x «^^«), (4) 

where in (4), H should correspond to the minimum value of the energy density 
as W{C) involves vacuum expectation value. 

The classical energy density for the background (2), in the Euclidean 
formulation, is £ — This energy density has a minimum S — sit H — 
and so W{C) in (4) does not give the area law. In order to realize the area 
law from (4), the minimum energy density should correspond to H ^ 0. 

Savvidy [1] has studied the quantum 1-loop effective energy density which 
has a minimum lower than the above classical minimum and for which H 0. 

However, Nielsen and Olesen [4] pointed out that the 1-loop effective energy 
density in the background (2) had an imaginary part, stemming from the 
lowest Landau level and so the vacuum (ground state) of such a model is 
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unstable. Various attempts were made to circumvent this sensitive issue 
which inhibited the progress. Instead of using the background (2), constant 
non-Abehan background — Q-^ Af — KSf was tried [5] and the said 
instabihty persisted. All these calculations were performed in the Gaussian 
(keeping only the terms quadratic in quantum fluctuations) approximation. 

We [6] have reexamined this important issue by retaining all the terms 
in the quantum fluctuations. Besides quadratic terms, there are terms cubic 
and quartic in quantum fluctuations. The detail of these calculations are 
given in Ref.6, in which the effective energy density has been shown to be 
real. Briefly, the Euclidean functional integral for SU{2) pure Yang-Mills 
theory 

S = J d'x{-\F;;^F;,}, 
F;, = d.At-d.Al + ge'^'^AlAt, (5) 
is expanded around the classical background A^ in (2) as 

Al = Al + a;, (6) 
and the quantum fluctuations are taken to satisfy the 'background gauge' 

Dfal = 0. (7) 

This gauge choice is important. First of all, there is no Gribov ambiguity 
[7] in using this background gauge. It has been shown by Amati and Rouet 
[8] that the multiphcity of classical solutions satisfying the gauge condition 
is an irrelevant issue for quantizing non-Abelian Yang-Mills theories in the 
background gauge and an unambiguous generating functional is now possible. 
The correct treatment of the zero modes of the 1-loop operator gives the 
background gauge relative to the classical solution. See also [9]. The crucial 
point is that under a gauge transformation U, the quantum fluctuations in 
(6) transform homogeneously, namely, = a^t°-,aij, — >■ Ua^U~^ [9]. Second, 
with the background gauge (7), and using (3), we have D^^(A|^ + a^) = and 
so the 'Maximal Abelian Gauge' or Abelian projection is realized for the full 
gauge field -|- a^. 
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Now using (6) and (7) in (5), the unambiguous Euclidean generating 
functional Z becomes, 

z = (8) 

with 

where 

Qac ^ (5f + 2^e«^^F;,. (10) 

In arriving at (9), we have introduced the gauge fixing and the Faddeev- 
Popov ghost Lagrangian for the background gauge (7) and integrated the 
ghost fields, resulting in the last term in (9). The expansion in (9) is exact. 
The purpose of writing 5" in the form above is to isolate the stable and 
unstable modes of 0^^. 

For the Sawidy background, ©44 = ©33 = D^^D^, so that their contri- 
butions to r cancel the ghost contribution. Further the non-vanishing ©'s 
are ©"J for i,j — 1, 2. Their eigenmodes and eigenvalues are: 

al±ial : ki^ + + k^^ + k^^ {plane waves), 

{al + ial) - i{al + ial) : (2n + l)gH + 2gH + kl + kj, (stable) 

{a\ - ial) + i{al - ial) : {2n + l)gH + 2gH + kl + kj, (stable) 

(a.l + ia.l) + i(al + ial) ■ (2n + I) gH - 2gH + kl + k^, (unstable) 

(a\ — ia\) — i(cL{ — icL^) : (2n + l)gH — 2gH + kl + kl, (unstable). 

The last two eigenvalues become negative when n — and for low mo- 
menta. These are called the 'unstable modes'. As we encounter logarithm of 

the eigenvalues, in the quadratic approximation, negative eigenvalues make it 
imaginary and hence the effective energy density becomes complex indicating 
vacuum instability. This in the Gaussian approximation. 
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The stable modes (the first two eigenvalues and the last two with n ^ 0) 
can be safely treated in the quadratic approximation. The contribution from 
the stable modes (see Ref.6 for details) are found to be 



967r2 



{£og{f)+C}, (11) 



where C is a real (infinite) constant and //^ is a dimensionful constant intro- 
duced to render the argument of the logarithm dimensionless. 

For the unstable modes, we [6] considered the full action in (9). The 
unstable modes involve the Lorentz indices 1 and 2 and the SU (2) indices 1 
and 2, because the classical background (2) is in the third color direction and 

so the cubic term in (9), namely, e'*'^'^(Z)"^apa^/i^ vanishes. The quartic term 
in (9) for the unstable modes is found to be |(|anP)^ where Gu is the unstable 
mode. The functional integral Z for the unstable modes is evaluated in Ref.6 
and from this the finite part of the unstable mode contribution to the energy 
density is found to be 

- "Sf %A (12) 

where 



dc' e-{^'^(^3 (13) 

The integral / is convergent irrespective of whether k'^^ + A;4^ is < or > 1. 
Further / is real, finite and independent of H. Adding (12) to (11) and 
including the classical energy density, the effective energy density is found to 
be 

where C includes the second term in (12) along with C in (11). The real 
constant C is then fixed by Coleman- Weinberg normalization ■^j2£\gH=ii? = 
I as — |. Thus the effective energy density for SU{2) Yang- Mills theory in 
Sawidy background becomes 
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This is real. The above result is non-Abelian gauge theory effect. 



Quarks (fermions) can be added by minimally coupling them with the 
background (2) and functionally integrating ip and ^ in Z. The only change 
is the replacement of 11 in (15) by (11 — Nf) for Nf quark flavors. The 
prefactor ^\g^2^ can be obtained from group theory considerations as well. 
Extending to SU{3), this factor becomes ^^g^^- 

In contrast to the classical energy density, the effective energy density in 
(15) has a minimum at non-zero H. Prom (15), we have 

dS _ 1 IV 

d^S 11 (/2 



> 0. (16) 



The energy density has a minimum. When quarks are included, in order to 
have a minimum energy density, Nf < 11 for SU{2) or Nf < 16 for SU{3). 
The minimum occurs when 

H = y e (17) 

with 11 appropriately replaced when quarks are included. Thus vacuum 
expectation value H ( corresponding to the minimum of S) is not zero which 
gives the Wilson loop the 'area law' and hence confinement.'^ 

The minimum energy density using (16) in (15) is 

£min - gg^2 ' (18) 



^In the case of QED, as photons do not have self-interactions, the electrons alone 
contribute to the effective energy density in the constant magnetic field background. In 
this case, the effective energy density will be £^ = ^ — A{e'^H'^){£og(^^^ — 5} where 
A is a positive constant. This energy density has a maximum at non-zero vahic of H. 
Interpretting the maximum of the energy density to correspond to excited state, the area- 
law gives confinement of electrons by a linear effective potential, in a constant magnetic 
field when the electrons are in the excited state, akin to magnetic confinement of plasma 
state. 
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which is lower than the classical minimum. This is the energy of the vacuum 
in the pure SU{2) Yang- Mills theory. 

The result that the minimum of the energy density occurs when H ^ 
(17) and F^2 = -^j imply that the vacuum expectation value (-^^2-^1^2) 7^ 0- 
This indicates that {g^F^i^FI^^) 7^ 0, the occurence of 'gluon condensate'. In 
this case 

(/C^;.) = 2///^,„ = 2^^'e-^. (19) 

Instead of using the strong coupling g which runs, we use the result from the 
Charmonium decay analysis, (^'^F^j.F^^) ~ 0. 5G eV^. Then, gH 0.5GeV^. 
With this estimate 

W{C) - "'■^'^ = e" (20) 

where a — ^ — 0.25GeV^. It is well known that the 'area law' corre- 
sponds to a linear potential and in the leading order V — aR where R is the 
separation of static quark and anti-quark.^ 

For a linear potential V — kr, the non-relativistic potential model calcu- 
lations give the cc bound states for k — 0.272GeV'^ which agrees with our 
estimate of a as 0.25GeV'^. 

Thus, the stable vacuum in the chromomagnetic background is very much 
indicative of confinement, giving in the leading order the linear potential 
whose parameter k is satisfactorily obtained. 

^The Wilson loop in (4) involves the area in the (a; — y) plane for the Savvidy back- 
ground (2). Following M. Baker, J. S. Ball, N.Brambilla, G.M. Prosper! and F.Zachariasen, 
Phys.Rev. D54 (1996) 2829, the closed loop is defined by the quark (antiquark) trajec- 
tories Zi{t){z2{t)) running from yi to x\ {x2 to ^2) as t varies from ti to tf. The quark 
(antiquark) trajectories are the world lines Ci (and C2) running from ti to tf (tf to ti). 
The world lines Ci and C2 along with two straight lines at fixed time connecting yi to y2 
and Xi to X2 then make up the contour. Parameterising the quark (antiquark) trajectories 
as: yi{-§,-'^)]Xi{-§,^);x2{§,^);y2{§,-^), where ti = -^,tf = ^ and R is the 
separation of the quark from antiquark, the area of the loop is RT. So, W{C) = e~*T"^-^. 
The potential is V = ^logW{C) in the limit T ^- 00 and so F = ^ii. 
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Now we consider the SU{2) Yang-Mills theory at finite temperature. In 
the studies of the Savvidy vacuum at finite temperature in the Gaussian 
Approximation, the effective energy density involved a temperature depen- 
dent imaginary part [10]. This inhibited the progress of examining the phase 
transition. We [11] have extended our zero-temperature studies (including 
the cubic and quartic terms) to finite temperature with chemical potential 
for the SU (2) gauge bosons. 

A chemical potential for massless non-Abelian bosons has been intro- 
duced in [12,13], by observing that there are conserved color charges = 
/ d^xj^; = /"-^"AIFI;^. For SU{2), one chooses Q^. The grand canonical 
partition function will now have iJ,Q^ in the hamiltonian. This leads to the 
result of using Aq — —ijiS"-^. This is not possible for Abehan gauge bosons. 

The role of chemical potential as a constant term in Aq is similar to 
the use of Polyakov loop specified by a constant A^ field in the third color 
direction. Now the Savvidy background becomes 

Ai = .0}, (21) 

which gives Ff^ — H and which solves the classical equation of motion. The 
background covariant derivative (Euclidean) now is 

Df = 8x5"^ + ge'^^Al^ ^Jie'^%x, (22) 

where v\ = (1, 0, 0, 0). 

Once again we have isolated the unstable modes and treated them includ- 
ing the cubic and quartic terms in the fluctuations. The involved calculation 
(see Rcf.ll for details) leads to the result that the effective energy density is 
real. The result for the effective energy density, including the zero tempera- 
ture contribution is 

+ Ml!g^{_|y,(^,^) + X,(^,^) 
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+ 2 ^ y2^H^Xi(v/2^H^^£^^//)}, 



(23) 



n=l 



where Yi, Ki are modified Bessel functions. Setting /3 = and ii = by/gH, 
the temperature dependence is written as 



TT 



1 ^ cos(a6^) , TT^^ , ^. . ^, 



{gH f 45a4 7r2a ^ £ '2 

oo 

+ 2 E V2n + lKi(a£V2n + 1)}. 

n=l 



(24) 



In [11], we have plotted (^p- with T in units of \/gH. The parameter 
h involves the chemical potential. For 6 = 0, zero chemical potential, the 
variation is smooth. At high temperatures, the behaviour is like that of non- 
interacting rclativistic gas. For h = 1,2,3 the variation shows a minimum and 
then smooth rise. A non-zero chemical potential thus triggers deconfinement 
phase transition. 

Deconfinement occurs for 6=1 around T j \J gH 0.4 and for b — 2 
around T/y/gH ~ 0.7. We have earlier identified the string tension o" as ^ 
and so our results give for the deconfining temperature ~ 0.5656 for 6 = 1 
and 0.9899 for b — 2. It is interesting to compare with the lattice studies. 
For SU{2), Lucini, Teper and Wenger [14] find the deconfining temperature 

~ 0.709; the agreement is satisfactory for 1 > 6 < 2. As the chemical 
potential fi = b\/gH, using gH = 0.5 GeV^, the lowest value for the chemical 
potential triggering deconfinement is 0.7 GeV < /i < 1.41 GeV. 



To summarize, the stable chromomagnetic vacuum for SU{2) Yang-Mills 
theory found in [6] gives a model for confinement using Wilson loop and hence 
a linear potential (in the leading order) for the quark-antiquark interaction. 
The coefficient k in this potential is ~ 0.25 GeV"^, in satisfactory agreement 
with non-relativistic potential model for charmonium. At finite temperature, 
the real effective energy density found in [11] is used to obtain estimates of 
the deconfining temperature and this reasonably agrees with the lattice study 
for SU{2). 
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